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We study the dynamical properties of spinless fermions on the checkerboard lattice. Our main 
interest is the limit of large nearest-neighbor repulsion V as compared with hopping |f|. The spectral 
functions show broad low-energy excitation which are due to the dynamics of fractionally charged 
excitations. Furthermore, it is shown that the fractional charges contribute to the electrical current 
density. 
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I. INTRODUCTION 

Excitations with fractional charge have fascinated 
physicists for some time. An early investigation by Su, 
Schrieffer and Heegeiii dates back to 1979 and deals 
with the chain molecule trans-polyacetylene (CH)n. This 
polymer has two ground states A and B with alternat- 
ing single and double bonds between C atoms. In one of 
them the double bond is on the right side of a given C 
atom while in the other ground state it is on the left side. 
A soliton or kink excitation has the property of an inter- 
face between the two ground states, i.e., on one side of 
the soliton we have alternating bonds like in ground state 
A and on the other side the bonds are arranged like in 
ground state B. Shifting the soliton along the chain does 
not change the energy. Su et al. showed that depend- 
ing on the occupational number of the kink state which 
can be either empty, singly or doubly occupied, the ex- 
citation has a spin or charge only. In other words, there 
exist excitations with spin-charge separation but no ex- 
citations with fractional charge yet. This changes when 
trans-polyacetylene is heavily doped. It turns out that at 
certain doping concentrations of either particles or holes 
there are not only excitations with charge and spin sep- 
arated but their charge is a fraction v of the electronic 
charge e onlyi^ In the simplest case v = ±1/3, ±2/3 etc. 
It should be mentioned that this phenomenon does not 
require interactions between electrons but occurs within 
a single- (or independent) particle description. It involves 
lattice degrees of freedom though, since a double C = C 
bond has a different length than a single C — C bond. 
Doped trans-polyacetylene serves as an example of a one- 
dimensional system with fractional charges. 

Fractional charges are also found in two dimensions 
in the much celebrated fractional quantum Hall effect 
(FQHE). It was Laughlin^ who introduced this concept 
here. He explained with it the behavior of electrons in 
high-quality GaAs/AlGaAs heterostructures in an ap- 
plied magnetic field B normal to the planes. When 
the field deviates slightly from B = nhc/ev, where n 
is the density of the electrons and i' is a rational frac- 



tion with odd denominator, excitations with charge ±j^e 
occur. As usual h is Planck's constant. The excitations 
with fractional charge ve are based on a new type of cor- 
related ground state called the Laughlin state. It has 
been also demonstrated that the excitations fulfill frac- 
tional statistics, i.e., when two of them are exchanged the 
phase changes by e^'^^ . In distinction to the previously 
considered trans-polyacetylene, electron-electron interac- 
tions are crucial here. In an appHed magnetic field, cor- 
relations become strong since the kinetic energy of the 
electrons is reduced. 

The question was left open whether or not fraction- 
ally charged excitations exist in two or three dimensional 
systems without magnetic fields. In Ref. jj| it was sug- 
gested that in a pyrochlore lattice, a prototype of a three 
dimensional structure with geometrical frustration, exci- 
tations with charge e/2 do exist. The pyrochlore struc- 
ture is found, for example, in the transition metal com- 
pound LiV2 04. In that compound, vanadium has a half- 
integer valency and electronic correlations are strong as 
indicated by a large 7 coefficient in the low temperature 
specific heat C — 7T. Although LiV2 04 has motivated 
work of the kind presented here, we do not claim that the 
theory of fractional charges applies to that specific ma- 
terial. Generally, prerequisites for fractional charges in 
the considered model are strong short-range correlations 
and certain band fillings. In order to study a system of 
this sort in more detail, a checkerboard lattice is consid- 
ered which also allows for numerical investigations. The 
checkerboard lattice can be thought of as a projection 
of the pyrochlore lattice onto a plane. Although there 
are differences in the physics of the two systems due to 
the different dimensions, one can learn much from the 
simpler two dimensional system. Numerous studies of 
the pyrochlore and checkerboard lattice were done and 
reported in Refs. |i5.6 7 8] for spin degrees of freedom 
and J3llflll3| for charge degrees of freedom. The ques- 
tion has not been posed yet how one would be able to 
detect these excitations with fractional charges in case 
they exist. It is the purpose of the present investiga- 
tion to shed some light on that topic by calculating spec- 




Figure 1: Checkerboard lattice (periodic boundary conditions 
are used for the calculations). Lattices with (a) 18 sites, (b) 
32 sites and (c) 50 sites are considered. 



tral functions and optical conductivity of fully spin polar- 
ized or, equivalently, spinless fermions on a checkerboard 
lattice. It is demonstrated that the fractionalization of 
charge leads to characteristic features which are absent 
when fractionalization is forbidden. Perhaps the present 
investigation may sharpen the attention of experimen- 
talists when doing photoemission experiments on, e.g., 
systems with spinel structure when the filling factors are 
right. Also the experimental progress in the generation 
of optical lattices should be noted. Recently it has been 
reported that 3D optical lattices can by generated and 
filled with either bosons or fermions to simulate relevant 
Hamiltonians, e.g., in Ref. [121 and citations therein. 

The paper is organized as follows. First we introduce 
in Section^ the model Hamiltonian that is used for our 
calculations and show how under certain conditions frac- 
tionally charged particles (fcp's) arise on a checkerboard 
lattice. In Section lTTTl we briefly dicuss the numerical de- 
tails and introduce the approximations we applied. The 
calculated spectral functions and optical conductivities 
are presented in Section ITvl Finally, the results and con- 
clusions are summarized in Section Fvl 



II. HAMILTONIAN AND LATTICES 



The model Hamiltonian 



H = -t 



E(4s+h.c.) + ^E 



n.inj 



(1) 



(^J> 



(i.j) 



describes spinless fermions on a checkerboard lattice (see 
Fig. with nearest-neighbor hopping amplitude t and 
nearest-neighbor repulsion V. The operators c^icj) anni- 
hilate (create) a particle on site i. The density operator is 
Hi = c\c^. To avoid boundary effects, which are expected 
to be strong, we choose periodic boundary conditions in 
actual calculations. 

The band structure of non-interacting {V = 0) parti- 
cles on the crisscrossed checkerboard lattice consists of a 




Figure 2: (color online) (a) Example of an allowed configura- 
tion of a checkerboard lattice at half filling with possible low- 
order hopping processes, (b)-(d) Adding one particle leads to 
two mobile defects, marked by triangles. Arrows mark mo- 
tions of particles. 



dispersive and a flat band 

e"(k) = -2t-4icos 
e+(k) = 2i, 
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where a is the lattice constant. In passing, we mention 
that in the presence of interactions correlations are gener- 
ally strong in the flat band regimes. The kinetic energy 
does not counteract optimization of repulsions in that 
regimei^ In the following we assume t > 0. Note that 
our coordinat e sy stem is rotated by 45° relative to that 
of, e.g., Ref. |lj|. The resulting difference in boundary 
conditions can lead to noticeable numerical differences in 
particular for very small cluster sizes. 

Our main interest is in the regime t/V ^ 1, at half flll- 
ing. When t — the ground state manifold is macroscop- 
ically degenerate: Every configuration with exactly two 
particles on each crisscrossed square is a ground state. 
This is the two dimensional equivalent of the tetrahedron 
rulai^ (or Anderson rule) . We refer to configurations ful- 
filling the tetrahedron rule as "allowed configurations". 
The macroscopic degeneracy is lifted for small but finite 
hopping amplitude t. 

In the following we discuss processes up to order t^/V^. 
They can be classified as self-energy contributions and 
ring exchanges (see Fig. El^a)). The self-energy contri- 
butions are identical for all allowed configurations and 
therefore the ground state degeneracy is not lifted by 
them. The total amplitude of ring exchange around 
empty squares is proportional to t'^ /V and vanishes for 
fermions. However, the macroscopic degeneracy is lifted 
by ring exchanges ~ t^ /V'^ around hexagons. Note that 
for fermions the signs of the corresponding matrix el- 
ements depend on the sequence in which the sites are 



enumerated jii^ 

Placing one additional particle with charge e onto an 
empty site leads to a violation of the tetrahedron rule on 
two adjacent crisscrossed squares (see Fig.[2Ib)). The en- 
ergy is increased by 4F since the added particle has four 
nearest-neighbors. Particles on a crisscrossed square with 
three particles can hop to another crisscrossed square 
without creating more violations of the tetrahedron rule, 
i.e., without increase of repulsive energy (see panels (c) 
and (d) in Fig.|2l. By these processes, two local defects 
(violating the tetrahedron rule) can separate over large 
distances thereby gaining additional kinetic energy of or- 
der t. This is observed in our calculations. We will see 
below in the context of Fig. that the bandwidth from 
an artificially restricted calculation in which the defects 
are not allowed to separate is much smaller than that 
from the unrestricted calculation. 

These findings motivate the introduction of the follow- 
ing effective Hamiltonian that acts only on the subspace 
of allowed configurations with a given number of viola- 
tions (e.g., two violations in Fig.|2l(b)-(d)): 

H,„ = .g^(-l)"oA^)(^| + |^)(^|+h.c.) 
- i E ^ (els + h.c.) V . 



(3) 



It includes the lowest order ring exchange process with 
g ~ t'^/V^ (compare B.ei. UM . We will refer to it as the 
t-g model. Note that the sign of the matrix elements de- 
pends on the configuration, i.e, for row wise enumeration 
of sites a process has a minus sign whenever the site in the 
center of the corresponding hexagon is occupied {uq — 1) 
and a plus sign when it is empty. The hopping t is pro- 
jected onto the manifold of allowed configurations plus 
fixed numbers of violations by the projector V. By intro- 
ducing the independent parameters t and g the effect of 
ring exchange onto the dynamics of fractional excitations 
can be estimated. The tendency for confinement can be 
studied even on rather small clusters by increasing the 
ring exchange strength g. 

One can expect that for an infinitely large system 
the two defects can propagate freely and independently 
throughout the system. In that case they should be con- 
sidered as two elementary excitations each having a dis- 
persion e(k) ~ t and each carrying the fractional charge 
e/2. Similarly, the repulsion-energy increase of AV is 
split between the two elementary excitations. We will 
refer to these excitations as fractionally charged parti- 
cles (fcp's). It should be emphasized that they are not 
quasiparticles in the sense of Landau's Fermi liquid the- 
ory because they are not adiabatically connected to inde- 
pendent fermions. Analogously, taking one particle out 
of the system generates two fractionally charged holes 
(fch's), i.e., two crisscrossed square with one particle only. 
Furthermore, quantum fiuctuations generate fcp and fch 
pairs. Starting from an allowed configuration, when one 



particle is hopping to an empty neighboring site it creates 
one crisscrossed square with three particles (fcp) and one 
crisscrossed square with one particle (fch) only. 

Whether or not two defects actually do separate com- 
pletely is a very subtle problem. The answer will depend 
on lattice type, dimensionality, spin vs. spinless parti- 
cles, etc. The answer may also differ for zero and finite 
temperatures T and depend on defect concentration. In 
any case, the key observation is that whenever two de- 
fects with charge ±e/2 separate they change the back- 
ground configuration. Infinite defect separation (decon- 
finement) is more plausible for liquidlike ground states 
than, for ground states with, e.g., broken translational in- 
variance. In the latter case, restoring forces occur which 
can lead to linear confinement, analogously to quark con- 
finement known from QCD. Thus a thorough analysis of 
the ground state of the undoped, i.e., half filled system 
is a vital step towards answering the above questions. 
First numerical results for an effective model on finite 
systems up to 64 sites have been presented in Ref. [l4J . 
A two-fold degenerate ground state with stripe order was 
found for the particular model on the checkerboard lat- 
tice. Nevertheless, various observations suggested that 
it would disappear in the thermodynamic limit. Recent 
calculations and analytical considerationsifi indicate that 
this supposition, i.e, that the order would disappear in 
the thermodynamic limit, was erroneous. The presence 
of an ordered ground state leads probably to weak con- 
finement of fcp's with a characteristic length of L lat- 
tice spacings. A simple estimate is given by Lg ~ t, i.e, 
L ~ t/g — V'^/t'^ lattice spacings. For the used parame- 
ters, this leads to a large characteristic length of several 
hundred lattice spacings (weak confinement). Note that 
the 3D pyrochlore lattice appears at present to have a 
deconfined phase in which the fcp's can separate^»Siil. 
The present work focuses on the dynamics of a weakly 
doped system (half filling plus one extra particle) by cal- 
culating the spectral functions and optical conductivity. 
We stress that the numerical simulations are limited to 
small clusters. As a consequence, as long as clusters of 
size L ^ t/g — V^/t'^ ^ 1 are out of reach for the exact 
diagonalization, we can not hope to find simple finite size 
scaling relations towards the true thermodynamic limit. 
The present study should be considered as a contribution 
towards the understanding of the local behavior of fcp's. 
In particular it is impossible to distinguish between a free 
or weakly bound pair of fcp's. 



III. NUMERICAL DETAILS 

Obviously, conventional approximation schemes such 
as mean-field theories and Green's function decoupHng 
schemes are not able to describe the strong local corre- 
lation expressed by the tetrahedron rule. Furthermore, 
until now neither a creation operator formalism nor a 
field theoretical description for the fcp's have been de- 
rived. This suggests numerical studies. Unfortunately, 



the fact that we are deahng with fermions rules out the 
use of standard Monte-Carlo techniques. Thus we have 
chosen the exact diagonalization of the Hamiltonian (QJ 
for small finite lattices, even though the numerical effort 
increases exponentially with system size. 

DiagonaHzation within the full Hilbert space was done 
for a lattice containing 18 sites. Larger systems can be 
considered if we restrict ourselves within the correspond- 
ing Hilbert spaces to certain low-energy sectors. They are 
defined by the number of violations of the tetrahedron 
rule. For the actual calculations we used the smallest 
possible Hilbert spaces allowing for the dynamical pro- 
cesses that we are interested in. These consists in the un- 
doped case, i.e, at half filling of the allowed configurations 
and those with one additional vacuum fiuctuation present 
(one fcp and one fch). In the doped case, i.e, with one 
particle added to the system, the configurations account 
for no other violations of the tetrahedron rule than those 
that are due to the added particle (two fcp's). We refer to 
the space spanned by these selected configurations as the 
"minimal Hilbert space" in each case. The dimensions are 
strongly reduced as compared with the full Hilbert space 
(see Table HJ. Calculations for the doped system have 



18 sites 



32 sites 



50 sites 



Half filling 


Full Hilbert space 


48 620 


601 080 390 


1.2641 10^"* 


Allowed states 


68 


2970 


67 832 


Minimal (1 fluct.) 


2 228 


168 858 


16 178 232 




System doped with one 


particle 






Full Hilbert space 


43 758 


565 722 720 


1.2155 10^'' 


Minimal (0 fluct.) 


1 323 


98 784 


8 698 450 


Confined (0 fluct.) 


612 


47 520 


1 695 800 


One extra fluctuation 


11 475 


2 435 808 


- 



polynomial expansionii^ We found essentially identical 
results for both algorithms. However, the implementa- 
tion of the Lanczos method turned out to be slightly 
faster. We rewrite first G{z) as 



where 



G{z) = {i;o\AA^\i^o){M 



A^^o) 



z- H 



\^o) 



./iMA^W) 



(5) 



(6) 



Then the state \(f)o) is taken as starting vector to gener- 
ate iteratively with the Lanczos algorithm an orthogonal 
basis for the Hamiltonian H . Using the tridiagonal form 
of the Hamiltonian with respect to the Lanczos basis and 
Kramers' rule, the expectation value can be easily rewrit- 
ten in terms of diagonal elements a„ and off-diagonal 
elements 5„ of the Hamiltonian in form of a continued 
fraction 



G{z) 



{MAA^i^^) 



z — ao 



bh 



(7) 



Well converged results were obtained already after several 
hundred iterations. 

There is, of course, always the question remaining in 
how far results for finite clusters are indications for the 
behavior of systems in the thermodynamic limit. Since 
we cannot go beyond 50 sites, the reader should be cau- 
tioned against over interpreting the results. 



IV. RESULTS 



A. Spectral functions 



Table I: Dimensions of the full Hilbert space and some relevant 
subspaces for different lattice sizes. 

also been done for two other subspaces. In one case, the 
two fcp's subspace is extended to three fcp and one fch 
in order to check the validity of the results for the two 
fcp subspace. We refer to this Hilbert space as "one extra 
fiuctuation". In the second case, we confine the two fcp's 
to adjacent criss-cross squares, i.e., particles do not split 
into fractional charges. Such calculations have the pur- 
pose to demonstrate that charge fractionahzation leads 
to qualitative different behavior. 

The spectral functions and optical conductivity of in- 
teracting many-particle systems are expectation values 
of the form 



G{z) = {iPo\A 



1 



z-H 



^t|^o) 



(4) 



and can therefore conveniently be calculated numerically 
by the Lanczos continued fraction methodic or kernel 



Direct insight into the dynamics of a many-body sys- 
tem is provided by the spectral function 



^(k,cj)=A-(k,cj)+A+(k,oj), 



(8) 



which is the probability for adding (+) or removing (- 
) a particle with momentum k and energy oj {h — 1) 
to the system. This function can be directly related to 
angular-resolved photoemission spectroscopy (ARPES). 
The particle contribution is defined by 



^+(k,cj) = lim 



(9) 



xIm(V'o |ck - 

u 

and the hole contribution by 



irj 



En-H 



41^0^) 



A-(k,w) 



= lim 

X Im«|4 



(10) 



1 



i-q-Eo + H 



c^.W,) 



Here {iI'q) is the ground state of the system with N par- 
ticles. A small value of 77 = O.lt is used for Lorentzian 
broadening. The operators c^ (cl) are obtained from the 
corresponding operators in real-space representation 



(11) 



where N^ denotes the number of k points in the extended 
BZ and the sum is taken over all lattice sites. The re- 
sulting integrated spectral density is 




D{uj) 



wX^^^^^ 



(12) 



For a system with a translationally invariant ground state 
or if an average over all degenerate ground states is con- 
sidered, D{uj) is conveniently calculated in real space rep- 
resentation as a local expectation value 

The integrated spectral density D{uj) for the Hamilto- 
nian 1^ is displayed in Fig. The different values of the 
nearest neighbor repulsion V, show the transition from 
an independent particle system to the strongly correlated 
limit. A rather small system is considered in order to al- 
low for a diagonalization in the full Hilbert space with 
reasonable computational effort. For ^ = the disper- 
sion is given by Eq. [2l D{u}) includes contributions from 
the allowed k-points (k = 27r/3 {n^jUy) with {nx,ny) e 
7?). The dispersive band e~ is completely filled and con- 
tributes only to D^{lo) while the flat band is empty and 
contributes exclusively to D^{lo). 

Adding a particle increases the energy by AV . Remov- 
ing it decreases the repulsion energy by 2V . Together 
this leads to an increasing separation of particle and hole 
part of the spectrum with increasing V and flnally to 
the formation of a gap. By this argument one expects a 
metal-insulator transition as a function of V . At small 
but flnite V the peaks are broadened as an incoherent 
background develops. But we are interested in the oppo- 
site limit of large V . 

In the limit of large V it is sufficient to calculate the 
spectral functions within the minimal Hilbert space, as is 
seen from the similarities with results for the full Hilbert 
space, shown for parameter values V = 16t and V — 32t 
in Fig. El In particular, the bandwidth as well as the 
gross structure remain unchanged. Therefore we may 
study much larger systems and compare results for dif- 
ferent lattice sizes. In particular, this enables us to ap- 
proach the question whether the two defects created by 
injecting one particle are closely bound to each other 
or not. Unfortunately, calculation on finite clusters can 
not distinguish truly between free particles and weakly 
bound pairs. For a first answer, we compare for differ- 
ent lattice sizes the results within the minimal Hilbert 
space with those from an artificially restricted calcula- 
tion keeping the two defects confined to two adjacent 
crisscrossed squares. From Fig. it is seen that without 
the restriction a broad low-energy continuum is obtained 
which is missing when the restriction is imposed. The 
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Figure 3: Integrated spectral density D{lj) of a vl8 x V 18 
checkerboard lattice for increasing nearest-neighbor repulsion 
V. The contributions of yl~(k, oj) is shown as shaded area. 
The insets show D{uj) calculated for rainiraal Hilbert spaces. 
Bandwidths and positions of the dominant features are almost 
unchanged. The peaks are broadened by choosing -q = O.lt. 



bandwidths are about 13t and 8t in the two cases. This 
suggests a simple interpretation. The dynamics of two 
separate fcp's having a bandwidth of « 6i each due to 
six nearest neighbors would yield the calculated lit while 
a confined added particle has a much smaller bandwidth, 
i.e., 8i. They are rather independent of the system size 
(e.g., ~ lit for the 18, 32 and 50 sites cluster in the free 
case) and therefore expected to extrapolate to the ther- 
modynamic limit. A similar observation holds true for 
A+(k,w) (see Fig. El where k = (0,7r/2)and k = (0,7r) 
are considered. Additionally one observes two interest- 
ing facts: (i) The states with large momentum k give the 
main contributions to the low-energy continuum in D{u). 
This is consistent with the picture that a particle with 
large momentum more efficiently lowers its kinetic energy 
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Figure 4: Integrated spectral density D{lu) for V — 25i calcu- 
lated for different lattice size. The left panel shows the data 
from "free" systems which include states with separated frac- 
tional charges. The right panel shows the data from systems 
where the two defect are "confined" and can not separate. The 
peaks are broadened by choosing r; = O.li. 



when decaying into fractional excitations than a particle 
with low momentum, (ii) >l+(k, w) shows in the "con- 
fined" calculation a peak-like feature at the low-energy 
edge that reminds one of the Landau-Fermi Hquid peak 
while nothing ahke is seen for the "free" case. 

A third observation deserves a comment. For vanishing 
momentum, the full spectral weight of A+ is contained 
in a single sharp i5-like peak near the center of the band. 
This suggests that |i/)^+^) — cj^=ol^o) is an exact eigen- 
state of the Hamiltonian H in the limit t/V -^ with 
if 1^^+^) = E\il;^~^^) as can be seen by evaluation of the 
following expression: 



H\^ 



N+l\ 



Hcl^oli^o) 



H,ci^o 



1^0 



Lo^lV'o 



H,cl^o iV'o) 



Eo\i^ 



N+l\ 



(13) 



The contributions of the kinetic and the repulsive energy 
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Figure 5: Spectral density A(k,aj) for V = 25i calculated for 
some fc-points. The upper panel shows data from the systems 
which include states with separated fractional charges. The 
lower panel shows the data from a system where the two defect 
are "confined" and can not separate. The peaks are broadened 
by choosing r\ — O.lt. 



to the commutator are given by 



-^ki„,Cp.^(, 



^£(k')4,Ci^,,cLo 

L k' 



and 



77 r^ 



V 

N 



E(4 



n, + c] n^ 



(14) 



(15) 



respectively. The latter is calculated by using the real 
space representation Ck=o = l/vjV^;c[. Since the 
ground state contains in the considered limit only config- 
urations that obey the tetrahedron rule, each empty site 
has exactly four occupied neighbors. Thus the sum over 
all nearest neighbors applied to the ground state leads to 









Y^.ti 



^^^Il^o) 



(16) 



Collecting everything yields 

H\^^+^) = (e(k = 0) + W + Eo) l^^+i). (17) 

The ground state wave function with one added particle 
of zero momentum is an eigenfunction of the Hamiltonian 
with energy £(k = 0)+4:V + Eo. Note that for the doped 
system no vacuum fluctuation is taken into account and 
thus one has to add a constant energy shift from the self- 
energy contributions (loss of ground state correlations) . 
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Figure 6: Integrated spectral density of a 32 sites cluster cal- 
culated for the effective t — g Hamiltonian for two different 
values of g. The peaks are broadened by choosing 77 — O.lt. 



Figure 7: Regular part of optical conductivity crreg(ii') of 
the vis X vis checkerboard lattice for increasing nearest- 
neighbor repulsion V. Left panels are showing the half filled 
case and the right half-filled plus one. The insets show the 
data for the same values within the minimal Hilbert space. 
The peaks are broadened by choosing rj = O.li. 



The integrated spectral density D{uj) of the effective 
Hamiltonian ^ shows a strong dependence on the ratio 
t/g. In the "physical" regime, corresponding to the pre- 
viously considered parameters with g ~ t^/V'^ = O.Olt, 
the integrated spectral density D{lj) (Fig. EJ shows qual- 
itatively the same features as the full Hamiltonian. If g 
is assumed to be equal to t, the broad continuum at the 
bottom of the spectral density vanishes and a sharp S- 
peak evolves instead. Note the similarities to D{uj) for 
the artificially confined situation. The spectral weight of 
the peak can be looked at as that of a Landau quasipar- 
ticle peak. This suggests the interpretation that the ring 
exchange term leads to charge order which is destroyed 
by the separation of the two fcp's. In the regime with 
g ^t {g = t^ /V'^) the diameter of the two bounded fcp's 
is larger than the considered system size and thus the ex- 
citations seem to be deconfined. An artificial increased g 
leads to much stronger confinement and the diameter of 
the bounded pair is small compared to the system size - 
leading to a finite weight of the quasiparticle peak. The 
huge spatial extend of the quasiparticle in the "physical" 
regime is expected to lead to very interesting effects, but 
this will be the subject of a separate investigation. 



B. Optical Conductivity 

The regular part of the optical conductivity cr,eg is de- 
fined by 
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Here |'0o') is the ground state of the system with N par- 
ticles and as before we use a value of 77 = 0.1 for broad- 
ening. The current operator is 



j^ = i[H,X] 
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(19) 



Calculations have been performed for the current den- 
sity along the x-axis, i.e. x = e^, in^l The regular part 
of the optical conductivity cr„g(w) is shown in Figured 
for increasing values of the nearest-neighbor repulsion V. 



For V = 0, one observes a peak-like structure. In the 
thermodynamic Hmit it should be positioned at w = 
and the shift to w ^ is a finite size effect. The same 
holds true for V=4t, which is expected to be still in the 
metallic phaseji^ With increasing V, the structures be- 
come broader. In the half-filled (undoped) case the com- 
plete weight is moved to larger uj and the Drude weight 
goes to zero. One expects a transition to an insulat- 
ing state. For large V the weight is distributed around 
u} = V. This corresponds to the energy that is needed to 
generate a fcp-fch pair which carry an electrical current. 

In the doped case one finds a different behavior. Fi- 
nite weight is found at a; = for arbitrarily large values 
of V and the charge is carried by two fcp's with charge 
e/2 each. A part of the weight is shifted to larger uj 
where additional fcp - fch pairs are generated as charge 
carriers. For V = 16t we compare the results based on 
the calculation in the full Hilbert space with those for 
the minimal Hilbert space. The features at small w are 
reproduced very well, but in the minimal Hilbert space 
vacuum fluctuations are absent and thus there is no cur- 
rent contributions from fcp-fch pairs. 

Now we want to investigate the influence of fractional 
charges on the optical conductivity (see Fig. [HJ. The 
main features are independent of lattice size. For the 
"free" system one flnds a bandwidth of nearly 13t which 
is the same as for the spectral functions. There are sharp 
peaks superposed onto a broad structure. The "conflned" 
system has a reduced bandwidth of nearly 8i and the 
broad structure is less pronounced. For the 32 sites clus- 
ter the contribution near w = is smaller than for the 
other two clusters. This is due to the particular shape 
of the clusters which lead to a lower degeneracy of the 
ground state than in the 32 sites case. 

In FigureElthe optical conductivity due to one moving 
fcp is shown. Here we keep one of the two fcp's flxed and 
let only the other one propagate. One observes a Drude 
peak and broad structure. Except for the Drude peak no 
other peaks are observed. The bandwidth is reduced to 
about half of the bandwidth of two fcp's. 



V. SUMMARY AND OUTLOOK 

We have studied numerically the dynamical properties 
of spinless fermions on checkerboard lattices and com- 
pared the results for different lattice sizes and models. 
For the full Hamiltonian a broad low-energy continuum 
is found and no Fermi Hquid peak is present, indicat- 
ing that an added particle decays into two fractionally 
charged excitations that separate over the whole flnite 
lattice. 

Considering the ring exchange g in the effective Hamil- 
tonian as an independent parameter (not flxed to t^/V^) 
allowed us to explore the regime where the fcp separa- 
tion is small compared to the system size. The spectral 
function for large g does not show a broad structure and 
develops instead a sharp peak. The existence of a quasi- 
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Figure 8: Regular part of optical conductivity (T,eg(tt;) for the 
doped system with V = 25t calculated for different lattice 
sizes. The left panel shows the data for "Free" systems which 
include states with separated fractional charges. The right 
panel shows the data for systems where the two defect are 
"confined" and can not separate. The peaks are broadened 
with rj = O.lf. 
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Figure 9: Regular part of optical conductivity cr,eg(a;) for the 
doped system with V = 25t. Only one of the fractional 
charged excitations is mobile while the other is fixed. The 
peaks are broadened by choosing tj = O.li. Arrows mark mo- 
tions of particles. 



particle peak shows that the added particle creates two 
bound fractionally charged particles with a small diame- 
ter. These findings suggest that for parameters V/t « 10 
quasiparticles with spatial extend over more than hun- 
dred lattice sites are formed. 

Finally the questions about the dynamics of fractional 
charges on different types of lattices remains. Consid- 
erable differences are expected for bipartite and non- 
bipartite lattices. Also the dimensionality of the lattice 
is crucial. For the 3D pyrochlore lattice one speculates 
about the existence of liquid phases where the fractional 
charges would be deconfined4iS*ii The informations we 
got from investigations of finite checkerboard lattices sug- 



gest that the possibly deconfined fcp's in the pyrochlore 
lattice would lead to broad features in spectral functions 
even for infinite systems. 
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